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We review various methods to investigate the statics and the dynamics of collective composition 
fluctuations in dense polymer mixtures within fluctuating-field approaches. The central idea 
of fluctuating-field theories is to rewrite the partition function of the interacting multi-chain 
systems in terms of integrals over auxiliary, often complex, fields, which are introduced by 
means of appropriate Hubbard-Stratonovich transformations. Thermodynamic averages like 
the average composition and the structure factor can be expressed exactly as averages of these 
fields. We discuss different analytical and numerical approaches to studying such a theory: The 
self-consistent field approach solves the integrals over the fluctuating fields in saddle-point 
approximation. Generalized random phase approximations allow to incorporate Gaussian 
fluctuations around the saddle point. Field theoretical polymer simulations are used to study 
the statistical mechanics of the full system with Complex Langevin or Monte Carlo methods. 
Unfortunately, they are hampered by the presence of a sign problem. In dense system, the latter 
can be avoided without losing essential physics by invoking a saddle point approximation for 
the complex field that couples to the total density. This leads to the external potential theory. 
We investigate the conditions under which this approximation is accurate. Finally, we discuss 
recent approaches to formulate realistic dynamical time evolution equations for such models. 
The methods are illustrated by two examples: A study of the fluctuation-induced formation 
of a polymeric microemulsion in a polymer-copolymer mixture, and a study of early-stage 
spinodal decomposition in a binary blend. 
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1 

Introduction 

Polymeric materials in daily life are generally multicomponent systems. Chemically 
different polymers are "alloyed" as to design a material which combines the favorable 
characteristics of the individual components [1]. Understanding the miscibility in the 
molten state is crucial for understanding and tailoring properties that are relevant 
for practical appUcations. Miscibility on a microscopic length scale is desirable, 
for instance, to increase the tensile strength of the material. Unfortunately, different 
polymers are typically not microscopically miscible, because the entropy of mixing 
is much smaller for polymers than for small molecules. The mixture separates into 
domains in which one of the components is enriched. The domains are separated 
by interfaces. Tuning the interface tension between the coexisting phases, and the 
morphology of the material on a mesoscopic length scale, is a key to tailoring 
material properties and has attracted abiding interest. 

One strategy to improve the miscibility and the interfacial properties has been 
to synthesize copolymers by chemically linking polymers of different type to each 
other. Added in small amounts to a homopolymer mixture, the copolymers adsorb 
at the interfaces, change their local structure, and reduce the interfacial tension [2, 
3, 4, 5J. Mixtures containing substantial amounts of copolymers form new types of 
phases, where different monomer types aggregate into mesoscopic "phase separated" 
domains (microphase separation). By varying the composition of the mixture and the 
architecture of the individual (co)polymers, one can create a wide range of different 
morphologies, corresponding to different materials, each with unique material properties. 

Unfortunately, copolymers are often rather expensive components. Another, 
cheaper, strategy for tuning the domain structure of blends takes advantage of 
the fact that in many practical applications, polymeric blends never reach thermal 
equilibrium on larger length scales. The morphology of the blend strongly depends 
on the kinetics of phase separation. Finer dispersed morphologies can be obtained by 
optimizing the processing conditions [6, 7]. 

One of the most powerful methods to assess such phenomena theoretically is 
the self-consistent field (SCF) theory. Originally introduced by Edwards [8] and 
later Helfand et al. [9], it has evolved into a versatile tool to describe the structure 
and thermodynamics of spatially inhomogeneous, dense polymer mixtures [10, 11, 
12, 13]. The SCF theory models a dense multi-component polymer mixture by an 
incompressible system of Gaussian chains with short-ranged binary interactions, and 
solves the statistical mechanics within the mean-field approximation. 

Much of the success of the SCF theory can be traced back to the extended size of 
the polymer molecules: 

First,the large size of the chain molecules imparts a rather universal behavior 
onto dense polymer mixtures, which can be characterized by only a small number 
of coarse-grained parameters, e. g., the end-to-end distance Rf. of the molecules, 
and the incompatibihty per chain where the Flory-Huggins parameter % [14] 
describes the repulsion between segments belonging to different components, and 
N is the number of segments per molecule. Rg sets the characteristic length scale 
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of spatial inhomogeneities. Since this length scale is much larger than the size of 
the repeat units along the backbone of the polymer chain, these systems can be 
described successfully by coarse-grained chain models. The parameters, R,, and 
X-/V, encode the chemical structure of the polymer on microscopic length scales. 
Indeed, comparisons between the predictions of the SCF theory and experiments or 
simulations have shown that the properties of many blends on large length scales 
depend on the atomistic structure only via the coarse-grained parameters Re and 
XN. The comparison works best if the parameters Re and x are considered as 
"black box" input parameters. Their determination from first principles or even from 
a well-defined microscopic model (e. g., atomistic force fields), is a formidable 
theoretical challenge: Even in the simplest case, the Flory-Huggins parameter x 
stems from small differences of dispersion forces between the different chemical 
segments. The value of the end-to-end distance Re of a polymer in a homogeneous 
melt results of a subtle screening of excluded volume interactions along the chain 
by the surrounding molecules, and Rp depends on the density and temperature. Thus 
X and Re cannot be calculated rigorously. If these coarse-grained parameters are 
determined independently (e. g., by experiments) and used as an input, the SCF 
theory is successful in making quantitative predictions. 

These arguments rationalize the success of a coarse-grained approach. The 
second reason for the success of a the SCF theory, which is a mean-field theory, is the 
fact that due to their extended shape, the molecules have many interaction partners. 
Let Q denote the polymer number density. The number of molecules inside the 
volume of a reference chain is then given by = gRf (where d denotes the spatial 
dimension). This quantity measures the degree of interdigitation of the molecules. In 
a dense three dimensional melt, Af is proportional to the number of segments N per 
chain, and it is refered to as the invariant degree of polymerization [15]. For systems 
which differ in J\f but are characterized by the same xN and Re the SCF theory will 
make identical predictions. The quantity Af plays an important role as it controls the 
strength of fluctuations. 

In this contribution, we shall discuss improvements which incorporate fluctuations 
into the SCF theory. Those fluctuations are particularly important (i) in the vicinity 
of phase transitions (e. g., the unmixing transition in a binary blend or the onset 
of microphase separation in block copolymers) or (ii) at interfaces, where the 
translational syrmnetry is broken and the local position of the interface can fluctuate. 
Both type of fluctuations can qualitatively change the behavior compared to the 
prediction of mean-field theory. The outline of the manuscript is as follows: First 
we introduce the Edwards model, which is employed as a coarse-grained description 
of a dense polymer melt throughout the article. Then, we briefly summarize the 
SCF theory for equilibrium properties and discuss various numerical strategies. 
Subsequently, we proceed to discuss possibihties to go beyond the mean-field 
approximation and obtain a description of the dynamics. Finally, we discuss some 
selected applications and close with a brief outlook on interesting problems. 
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The Edwards model for polymer blends 

2.1 

Gaussian model 

In dense binary blends, many of the interesting phenomena occur on length scales 
much larger than the size of a monomeric unit. Hence, a theoretical description can 
use a coarse-grained model, which only incorporates the relevant ingredients that 
are necessary to bring about the observed universal behavior: chain connectivity and 
thermal interaction between unlike monomeric units. The universal properties on 
length scales much larger than the size of a monomeric unit do not depend on the 
details of these interactions. One can therefore choose a convenient mathematical 
model, the Edwards Hamiltonian [8], in which the local microscopic structure enters 
only via three phenomenological parameters, Rg, x-^, and Af. 

In a dense melt, the excluded volume of the monomeric units is screened and 
chains adopt Gaussian conformations on large length scales. In the following, we 
shall describe the conformations of a polymer as space curves r(T), where the 
contour parameter r runs from to 1. The probability distribution Vlr] of such a 
path r(r) is given by the Wiener measure 



This describes the Gaussian chain conformations of a non-interacting chain. It is 
characterized by a single length scale. Re the end-to-end distance of the polymer 
chain. The structure is Gaussian on all length scales, i. e., the model ignores 
the structure on short length scales [16]: Self-avoiding walk statistics inside the 
excluded volume blob, rod-Uke behavior on the length scale of the persistence 
length, and the details of the monomeric building units at even shorter length scales. 
All this structure on the microscopic length scale enters the Gaussian chain model 
only through the single parameter Rg. This parameter (and its dependence on the 
thermodynamic state, i. e., temperature, density or pressure) cannot be predicted 
by the theory but is used as input. The Gaussian model provides an appropriate 
description if the smallest length scale Zmin of interest in the calculations is much 
larger than microscopic length scales. The coarse-graining length scale Imin should 
be at least on the order of a few persistence lengths. If this condition is violated, 
the Gaussian chain model can produce quaUtatively incorrect predictions (e. g., at 
low temperatures where both the persistence length length becomes large and the 
interface width becomes small). In this case, other chain models have to be employed 
that take due account of the local structure (e. g., the worm- like chain model [17], 
which captures the crossover from Gaussian conformations at large length scales 
to rod-like behavior on the length scale of the persistence length, or enumeration 
schemes [18, 19, 20, 21] which can deal with arbitrary chain architecture). 

In the following we consider a polymer blend of homopolymers of species 
A and B and/or copolymers containing both types of monomers. For simplicity. 




(1) 
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we assume that segments of type A and B are perfectly symmetrical and that 
all polymers have the same chain characteristics, in particular, the same chain 
length. The generalization to more asymmetrical situations is straightforward. We 
characterize the distribution of segments A and B along a chain of given type 
J (corresponding to a ^-homopolymer, a S-homopolymer or an A : B block 
copolymer) by the segment occupation functions 7ja(t), 7jb(t) [22], which take 
the value or 1 and fulfill 

IJA (t) + 7jb (t) = 1 for allJ and r . (2) 

Combining these with the chain conformations r(T), we can define a microscopic 
density (/)a(r) for segments of type a = Aoi B 



Here the sum ij runs over all nj polymers of type J, and rjj(T) denotes the 
conformation of the i^j^ polymer. The segment density is normalized by the polymer 
number density, g = J2j '"^j/^- 

In addition to the chain connectivity, the coarse-grained model has to capture the 
interactions between segments. In general, a compressible binary polymer mixture [23] 
exhibits both liquid-liquid immiscibiUty as well as liquid- vapor type phase separation 
(cf. Poly mer+ solvent systems: Phase diagrams, interface free energies, andnucleation 
in this issue). In the following, however, we regard only liquid-liquid unmixing into 
A-rich and i?-rich phases or microphase separation into domains comparable in size 
to the molecular extension. 

To describe the interaction in the polymer liquid, the interaction potential between 
segments can be decomposed in a short-ranged repulsive part and a longer-ranged 
attractive part (e. g. , using the Barker-Henderson scheme [24] or the Weeks-Chandler- 
Anderson decomposition [25]). The short-ranged repulsive contribution determines 
the packing and structure of the fluid. The length scale of fluctuations in the total 
density is set by the excluded volume correlation length ^ev, which also sets the 
length scale below which the chain conformations do not obey Gaussian but rather 
self-avoiding walk statistics. In a dense melt, ^ev is a microscopic length scale. On 
length scales larger than ^ev density fluctuations have already decayed, i. e., the fluid 
is incompressible. For ^cv ^ Umiy, we can represent the effect of the short-ranged 
repulsive interactions by an incompressibility constraint: 



where we have assumed that both segments occupy identical volumes. Of course, 
this incompressibility constraint can only be enforced after the coarse-graining 
procedure, i. e. , on length scales larger than • The short-ranged repulsive interactions, 
hence, do not enter the Edwards Hamiltonian explicitely, but only via the number 
density of polymers g. The dependence of g on the thermodynamic state {i. e., the 




(3) 



(I)a{y)+(I)b{y) = 1 



(4) 
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temperature, pressure, and composition of the mixture) - the equation of state - is 
not predicted by this computational scheme. 

The longer-ranged attractive interactions between the segments do not strongly 
influence the local fluid structure, but drive the phase separation. Again we assume 
that the range of the attractive interactions is small compared to smallest length scale 
^min of interest in our calculations. Then, the detailed shape of the longer-ranged 
potential does not matter either, and we can represent the interactions as zero-ranged. 
By virtue of the incompressibility constraint, there is only one independent bilinear 
form of densities and so we use as interaction free energy the expression 



kBTQxN / d'^r <t>A{v)ct>B{v) 



(5) 



The combination of Flory-Huggins parameter and chain length, xTV, parameterizes 
all the subtle interplay between liquid-like structure of the polymeric liquid and the 
small differences in the attractive interactions between segments of different types. 
It is worth noting that the Flory-Huggins parameter for neutral polymers without 
specific interactions is on the order 10"** — 10~^ [26], while the strength of the 
attractive interactions in the fluid are of order 1 fc^T. Hence the x parameter, which 
describes the difference in attractive interactions between the species, arises from a 
strong cancellation. Small changes in the thermodynamic state might perturb this 
balance. As a result, the x-parameter depends on the temperature, pressure and 
composition, as is often observed in experiments. 

Within the coarse-grained model the canonical partition function of a polymer 
mixture with nj chains of type J takes the form: 



1 



Y{jnj\ 



Y{jnj\ 



[]p[ri]P[ri]exp 



n n '^[^^^n^i^^^v 



S{(pA + '/>B - 1) 



MxN 



exp 



MxN 



(5('/>A + 0B-l) 



(6) 



The product runs over all n = ^ nj chains, and the functional integral 
X'[rj] sums over all possible conformations of the chain i. Within the Edwards 
Hamiltonian, the system is thus described by the following parameters: describes 
the incompatibility of the different polymer species, i?,. sets the size of the molecules, 
M is the invariant chain length, which describes the degree of interdigitation of the 
molecules, = nj/QV is the average density of type J, and 7jq, (t) the distribution 
of segments along the molecules of type J. The quantity M sets the scale of the free 
energy, but does not influence the SCF solution {i. e., the location of the extrema of 
the free energy, see Sec. 3). 
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2.2 

Hubbard Stratonovich transformation and fluctuating fields 

The partition function (6) describes a system of mutually interacting chains. Introducing 
auxiliary fields U and IT' via a Hubbard-Stratonovich transformation, one can 
decouple the interaction between the chains and rewrite the Hamiltonian in terms 
of independent chains in fluctuating fields. Then, one can integrate over the chain 
conformations and obtain a Hamiltonian which only depends on the auxiliary fields. 
Thermodynamic averages like density or structure factors can be expressed as 
averages over the field U and W without approximation. 

In the following we introduce a minimal transformation to decouple the interactions 
between polymer chains. Two different schemes have to be employed for the thermal 
interactions between the monomers and the incompressibiUty constraint. For the 
thermal interactions which give rise to the term {(pA — <^s)^ in the Hamiltonian, 
we use the Hubbard-Stratonovich formula 



exp[ax 



1 



+00 



dy exp 



xy 



(7) 



at each point in space and identify X = ^a — ^b, y = anda = vMxN /AR^. 

This introduces a real auxiliary field W. 

To rewrite the incompressibiUty constraint, we use the Fourier representation of 
the 5-function: 



1 /■+°° 

^(a; - 1) = — j dy exp [-iy{x - 1)] 



(8) 



using X = (j)A + (j^B and y = gU/2. This introduces a real auxiliary field U. 

The field W, which couples to the composition of the mixture, gives rise to a real 
contribution in the exponent, while the field U, which enforces incompressibility, 
yields a complex contribution. Replacing the incompressibiUty constraint by a finite 
compressibility 



+ - 1)' 



(9) 



where k denotes the isothermal compressibiUty of the polymer liquid, would also 
introduce a complex term in the exponential via the Hubbard-Stratonovich formula 
(cf. Sec. 4.4.2). 

Using these two expressions we can exactly rewrite the partition function as 

qVxN' 



1 



exp 

VUVWeyi^ 



Wv[v,]V[Yi\VU[v\DW[v] exp 



/ 



' X^XN 
Q[U,W\ 
kBT 



4 

iU 



(10) 
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which defines a Hamiltonian g[U, W] [27, 28] for the fields U and W, 



g[u,w] 



NkBT{V/Ri) 



VQ[Wa,Wb] , 1 



nj 



-/d^r 



4xiV 



(11) 



Here Qj \Wa , Wb] denotes the partition function of a single noninteracting Gaussian 
chain of type J 



Qj[Wa,Wb] = - 



'D[r]V[r] exp 
2?[r]P[r] exp 



■^0 a=A,B 

i J d^r ^iy„(r)^yjr) 



in the fields 



Wa = 



iU + W 



Wb = 



iU-W 



(12) 
(13) 

(14) 



Here we have normalized Q j by the volume V in order to make the trivial contribution 
of the translational entropy explicit. The dimensionless, microscopic single chain 
density is defined as 



yjr)=V / dT<5(r-r(r)) 



(15) 



where a = A or B describes the segment species and J denotes the type of 
polymer. In Eq. (11), we have omitted a term — (1/2V) / d'^r {iU — /2) . This 
is legitimated by the fact that g\U, W\ is invariant under adding a spatially constant 
field AU to U. Therefore, we can fix 



(16) 



To calculate the single chain partition function Qj, it is useful to define the end 
segment distribution qj{v, t), which describes the probability of finding the end of 
a chain fragment containing the segments [0 : t] and exposed to the fields Wa {ct = 
^, B) at position r [9]: 

qj{Y,t) = J V[r{t)]V[r{t)]6{r{t) - r) e' drj:^w^ir{r)hjM) (^^^ 
Similarly, we define 

ql{r,t) = J V[r{t)]r[r{t)]d{r{t) - r) e"/*' -^^S" ^"(■•M)^^"^. (18) 



The end segment distributions obey the following diffusion equation [9] 
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^ = fvV(M) 
-%^ = fvM(M) 



J2 WaJja{t)qj{r,t) 

a=A,B 

J2 W^a7Ja(i)3j(r,i) 



(19) 
(20) 



=A.B 



with the boundary condition (j'j(r,0) = 1 and (/^(r, 1) = 1, e. , the beginning of the 
chain fraction is uniformly distributed. For homopolymers the two propagators are 
related via qj{r, t) = qj{r, 1 — t). The solution of these equations yields the single 
chain partition function: 

(^J = ^j qj{r, 1) = ^ / d'r g^(r,0) = ^ / d'r qj{r,t)ql{r,t) Vi 

(21) 

We note that the fluctuations in the spatially homogeneous component W = 
(1/V) / d'^r W of the field W are strictiy Gaussian. To show this, we decompose 
W = W + W' and obtain: 



G[U,W] 



g\u,w'] 



MkBT{V/Ri) 



+ 



MkBT{V/Ri) 
^ [w + xN{^a-4>b)Y -^{4>A 



Ib? (22) 



The only effect of a spatially homogeneous field acting on ^-segments or B- 
segments is to introduce additional contributions to the chemical potentials jij - 
quantities which are immaterial in the canonical ensemble. 

By the Hubbard-Stratonovich transformation we have rewritten the partition 
function of the interacting multi-chain systems in terms of non-interacting chains 
in complex fluctuating fields ill + W and iU — W. In field theoretical polymer 
simulations, one samples the fields U and W via computer simulation using the 
above Hamiltonian (cf. Sec. 4.4). 

To calculate thermal averages of (f)A{^) — (f>B{j') we introduce a local exchange 
potential Afi that couples to the local composition 4>a — (f>B [28]: 

Z[Atx] ~ j^^^ f Y[v[vi]V[ri]VU[r]VW[v] 



Ujnjl 
exp 



/ 



VUVW exp 



W 



g[u,w,An] 



(23) 



The free energy functional g[U, W, Zi/i] takes the form 
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g[U,W,Afx] ^- , VQJ{^{^U + W + A^i),^{lU-W-An)] 



VMkBT{y/red) j nj 



gp.w] , 1 



_ , ^ ^^^^ , d'^r {An^ -2 A^iW) (24) 

where we have changed the dummy field W to W+Aji in the last step. Differentiating 
Eq.(24) with respect to Ajiijc) we obtain: 



, . - , 2 1 5Z\Aii] 

i(r)-0B(r))- ^ 



Q Z[Atj\ 5Aii{v) 



= -^m (25) 

/\([i=0 



^2 J[Z\/i] 5Aii{Y)5Aii{r') 



This expression relates the thermodynamic average of the composition difference to 
the average of the real field W . Similarly, one can generate higher moments of the 

composition: 

([0^(r) - <^B(r)][0^(r') - ^b{v')]) (26) 
4 1 5'^Z[An] 2Ri6{v - r') {W{Y)W{r')) 

The two expressions (25) and (26) allow us to calculate the physically important 
thermodynamic average of the microscopic densities and their fluctuations from the 
thermodynamics average and fluctuations of the field W. Although the Hamiltonian 
Q\U, W] is complex, the thermodynamic average of the microscopic densities and 
their fluctuations can be expressed in terms of the real field W [28]. 

Alternatively, one can calculate the average of the microscopic densities from 
Eq.(24) and obtains 

(<^^(r) - 0B(r)) = {c^\[Wa.Wb] ~ <P*bWa.Wb]), (27) 

where the 0* {a = A, B) are functionals of Wa = (iU + W)/2 and Wb = {iU - 
W)I2 (cf.Eq. (14)) given by 

- /P[r]P[r]e-v/d''-S.^a(rfe(--)v/dT^(r-r(T))7j^(T)) 
^Y.^j{4>'U^))- (29) 

,7 

This equation identifies 0* as the density of a-segments created by a single chain 
in the external fields Wa and Wb, averaged over all polymer types J. (We note 
that the average normalized microscopic density of a single chain in {(j)^'') in a 
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volume V is independent from V, cf. Eq. (15).) The thermodynamic average of 
the composition is simply the Boltzmann average of the corresponding single chain 
properties averaged over the fluctuating fields. The functional derivatives 4>'\{r) and 
(?!)|j(r) can be calculated using the end-segment distribution functions, qj{r, t) and 
qj{r,t), according to: 



dt qj{r,t) q\{v,t) Jjait) 



(30) 



While the thermodynamic average of the microscopic composition (j>a equals the 
average of the functional (/)* over the fluctuating fields Wa and Wb, or, alternatively, 
U and W, such a simple relation does not hold true for the fluctuations. Taking the 
second derivative with respect to Afi, we obtain: 

([^^(r) - Mr)][Mr') - M^')]) = mi^) - <P*B{^m*A{r') - 'A^(r')]) 
Rt I S^Air) SrBjr) ScP\{r) J.^^(r) 
^ \ 5WA{r') 6WB{r') ^W^s(r') (5H^^(r') , 



(31) 



The fluctuations of the physically relevant microscopic monomer densities (I)a and 
(j)B and the fluctuations of the Boltzmann averaged single chain properties </)^ and 
0^ due to fluctuations of the fields U and W are not identical. The additional term 
accounts for the single chain correlations [29]: 



SWpir') 



V 



{r)$%{r')) - {4>7Jr)){$f0{r'))) . (32) 



We emphasize that the additional term due to single chain correlations is of the 
same order of magnitude than the fluctuations of the physically relevant microscopic 
monomer densities itself. We shall demonstrate this expUcitly in the framework of 
the Random-Phase Approximation (RPA) in Sec. 4.2. 

To calculate thermal averages of ^!>A(r) -|- (^s (r) we introduce a spatially varying 
total chemical potential 5fi 

Z[5ii] ^ :pj-i^ / l[V[ri]V[vi]VU[r]VW[r] 



exp 



VUVW exp 



l4xJ 

g[u,w,s^i] 



W .2 



iU + 6fj. 



t>A + 4>B 



(33) 



Similar to Eqs.(24) and (24) we can rewrite the free energy functional Q\U, W, 
as: 

g[u,w,5^l] g[u,w] 



VAfkBT{V/Ri) VWkBT{V/Ri) 

VQj 



nj 



Vj 4xN 



(34) 
(35) 
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The first derivative wrt (5/i yields the average of the local monomer density 

2 1 6Z[Sfi] 



e Z[5n] S6ii{r) 



{4>*A + 4>*b) 



Sn=0 



1 



The second derivative yields the density fluctuations: 

4 1 S'^Z[5^i] 



J[\mA{v')^ ^Wb{v') ' SWB{r') ' 5WA{r') 



1 



(36) 
(37) 



(38) 
(39) 



Note that the last contribution in the equation above is similar to the single chain 
correlations in Eq. (31). The incompressibiUty constraintis enforced on the microscopic 
density (/jyi + . At this stage, (f>*y^ and (pg are only auxiliary functionals of the fields 
U, W proportional to the density distribution of a single chain in the external fields. 



3 

External potential (EP) theory and self-consistent field 
(SCF) theory 

The reformulation of the partition function in terms of single chains in the fluctuating, 
complex fields Wa — {iU + W) /2 and Wb = {iU — W)/2 is exact. The numerical 
evaluation of the functional integral over the fields U and W is, however, difficult 
(cf. Sec. 4.4), and various approximations have been devised. 

The functional integral over the fluctuating field U , conjugated to the total 
density, can be approximated by the "saddle point": The integrand evaluated at that 
function C/*[W^] which minimizes the free energy function G\U, W]. Carrying out 
this saddle point integration in the field U, we obtain a free energy functional 7i[VK]. 
In the following we denote this external potential (EP) theory , following Maurits 
and Fraaije who derived the saddle point equations heuristically [30]. This scheme 
still retains the important fluctuations in the field W, conjugated to the composition. 
Using an additional saddle point approximation for W, we neglect fluctuations in the 
composition and we arrive at the self-consistent field theory. 



3.1 

External potential (EP) theory: saddle point integration in U 

The fluctuations described by the two fields U and W are qualitatively different. This 
is already apparent from the fact that one, U, gives rise to a complex contribution to 
the field that acts on a chain, while the other, W, corresponds to a real one. The 
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field U couples to the total density ipA + (pB and has been introduced to decouple 
the incompressibility constraint. Qualitatively, it controls the fluctuations of the total 
density, and it does not directly influence the expectation viilue or the fluctuations of 
the composition (j)A — (pB- 

The saddle point value U* is given by the condition: 

5G[U\W] 

=0 (pA + (PB = ^- (40) 

Instead of enforcing the incompressibiUty constraint on each naicroscopic conformations, 
we thus only require that the single chain averages in the external field obey the 
constraint. We recall that the 0* are functional of Wa = {iU + W)/2 and 
Wb = {iU - W)/2 (cf. Eq. (28)), thus Eq. (40) impUcitly defines a functional 
U* \W]. Substituting the saddle point value into the free energy functional (11), we 
obtain an approximate partition function 



/ 



•Zep ~ / T^W exp 



n[w] 



(41) 



with the free energy functional 

n[w] Q[u\w] 



.-^l,,.YQi,U,',f- ,42) 

nj V J 4X-/V 

At this point, one has two possible choices for calculating the thermodynamic 

averages of the microscopic composition [28]: 

1. Although one performs a saddle point approximation in U, one can use the 
expressions (25) and (26), which relate the average of the field W and its 
fluctuations to the average and the fluctuations of the microscopic composition. 
The crucial difference is, however, that the fluctuations of the field W are now 
governed by the Hamiltonian H[W] of the EP theory and not by the exact 
Hamiltonian ^[?7, VK]. If the saddle point approximation is reliable, tj[?7, W^] 
will be well described by a parabola inU — U* [W], the fluctuations in U will be 
nearly Gaussian and will hardly influence the fluctuations of W. Therefore, the 
fluctuations of W with respect to the Hamiltonian 7i[II^] of the EP theory closely 
mimic the fluctuations of the field W with respect to the exact Hamiltonian 
g[U,W], i. e., (W^(r)T4^(r'))EP « {W{r)W{r')). By the same token, using 
Eq. (26), we ensure that the composition fluctuations in the EP theory closely 
resemble the exact thermodynamic average of fluctuations of the microscopic 
composition. 

2. One can introduce a spatially varying exchange potential dfj, and calculate the 
thermodynamic averages with respect to the free energy functional H[W] of the 
EP theory. This leads to expressions sinnilar to (27) and (31): 
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(0A(r)-0B(r))EP = (<^i-</'|j)EP (43) 
{[Mr) - Mr)][Mr') - <^s(r')])EP (44) 

Ri I <50i(r) , d<l>*^(y) 6^^) S<l>%{r) 



V:^\SWAir') SWeir') SWBir') ,5T.y^(r') / gp 

We refer to the last expression as the Uteral composition fluctuations in the 
EP theory. Likewise, one can use the functional derivative with respect to a 
spatially varying chemical potential Z\/z, which couples to 4>a + 4'b, to calculate 
the literal average of the total density. 

(<^A(r)+.^s(r))EP = (<^:i + </>^)EP = l (45) 

For the literal fluctuations of the total density in the EP theory (cf. Eq.(38)) one 
obtains: 

([<^^(r) + <^B(r)][0^(r') + .^B(r')])EP (46) 
= (['/'A(r) + + <f>W)])EP 

V , , " 

=1 =1 



This demonstrates that the saddle point approximation in \J enforces the incompressibiUty 
constraint only on average, but the literal fluctuations of the total density in the 
EP theory do not vanish. 

Unfortunately, the two methods for calculating composition fluctuations yield different 

results. Of course, after the saddle point approximation for the field \J is is not a 
priory obvious to which extent the EP theory can correctly describe composition 
fluctuations. In the following we shall employ Eqs. (25) and (26) to calculate the 
thermal average of the composition and its fluctuations in the EP theory. As we have 
argued above, these expressions will converge to the exact result if the fluctuations 
in U become Gaussian. 



3.2 

Self-consistent field (SCF) theory: saddle point integration in U and W 

In the self-consistent field (SCF) theory one additionally approximates the functional 
integral over the field W by its saddle point value. Using the functionals 



blWA,WB\ = - 



(cf. Eq. (28)), the derivative of the free energy functional 'H\yV\ (Eq. (42)) can be 
written as 
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6n[W] _ VUkBT W + xN{4>*A - <^b) 



5W{v:) 



2xN 



(47) 



where we have used the convention (16). Therefore tiie saddle point condition simply 
reads: 



(48) 



where (f)^ and (j)^ are functionals of the fields U and W. 

The two equations (40) and (48) define a self-consistent set of equations for the 
two fields W* and U* . Knowing the fields W* and U*, one can calculate the thermal 
average of the microscopic composition via: 



W* 

{4>A - (I)b)sCF = 



(49) 



3.3 

Alternative derivation of the SCF theory 

The SCF equations are often derived in a slightly different way by additionally 
introducing auxiliary fields also for the microscopic densities. For completeness and 
further reference, we shall present this derivation here. 

The starting point is the partition function, Eq. (6). Then, one proceeds to 
introduce collective densities <Pa, ^b, fields Oa, ^b, and a pressure field U to 
rewrite the partition function in the form: 



5{<Pa +<Pb- 1)S{<Pa 



rr—i I Tmri]V[vi]VM^B'Df2A'Df2B'DU 



exp 



exp 



/ 



g J d^r {+if2A - ^a) + if^B - ^b) } 
:F[$A,^B,!^A,flB,U] 



'D^A'D^B'Df2A'DnB'DU exp 



where we have used the Fourier representation of the (5-functions. The free energy 
functional T[<Pai'^b^ ^Ai ^b,U] is given by: 



TI^a.^b.^a^Qb^U] 



VAfkBT{V/Rd) 



VQj[tf2A,if2B] 
nj 



iU 



XN^A^B + — (^A + - 1) - i^A^A - if^B^B 



(51) 
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The saddle point approximation of the partition functions yield the SCF theory. Let 
us first make the saddle point approximation in the fields. Like before, we denote the 
values of the fields at the saddle point by an asterisk. 



J 

E 

J 



= = - > 0j = 



Q^'^A+^B = l (52) 

Note that the saddle point values, i[2\ and iQ*^ are real. These equations correspond 
to Eq. (40). Inserting these saddle point values into the free energy functional, we 
obtain a functional that depends only on the composition <?^: 

T*[^a] _ T[^A,l-^A,fi*A,0*B,U*] 

^/MkBT{V/Rf} Vj^kBT{V/Rt) 
^_Y7^^}W^^aMi}!Ib\M 

+ ^ / d'^r {xN^a{1 - ^a) - i^i*AWA]^A - ^^^bI^aKI - ^a)} (53) 

This functional can be used to investigate the dynamics of collective composition 
fluctuations (cf. Sec. 5). If we proceed to make a saddle point approximation for the 
collective density ^a, we arrive at 

:^lM=o^i(/2* -/2J)+XiV(^i-rB)=0 (54) 

which is equivalent to Eq. (48) when we identify W* =in^ — iQ*^. 



3.4 

Numerical techniques 

Both in the SCF theory and in the EP theory, the central numerical task is to solve 
a closed set of self-consistent equations one or many times. This is usually achieved 
according to the general scheme: 

(0) Make an initial guess for the fields U or {W, U). 

(1) From the fields, calculate the propagators q and (Eqs. (19 and (20)). 

(2) From the propagators, calculate a new set of fields. 

(3) Mix the old fields and the new fields according to some appropriate prescription. 

(4) Go back to (1). Continue this until the fields have converged to some value within 
a pre-defined accuracy. 
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The numerical challenge is thus twofold. First, one needs an efficient iteration 
prescription. Second, one needs a fast method to solve the diffusion equation. 

We will begin with presenting some possible iteration prescriptions. Our list 
is far from complete. We will not derive the methods, nor discuss the numerical 
background and the numerical stability. For more information, the reader is referred 
to the references. 

We consider a general nonhnear set of equations of the form 

g(x) = X. (55) 

The value of x after the nth iteration step is denoted x„. 

One of the most established methods to solve nonlinear systems of equation, 
which has often been used to solve SCF equations, is the Newton-Raphson method [31]. 
In this scheme, one calculates not only g(x„) from a given x„, but also the whole 
Jacobian matrix G„ = The (n + l)th iteration is then given by 

x„+i = x„ + [1 - G„]"^(g(x„) - x„). (56) 

If this method converges at all, i. e. , if the initial guess is not too far from the solution, 
it usually requires very few iteration steps. Unfortunately, it has the drawback that 
the Jacobian must be evaluated in every iteration step. If one has m independent 
parameters (fields), one must solve 2m diffusion equations. Therefore, the method 
becomes inefficient if the number of degrees of freedom m, i. e. , the dimension of x, 
is very large. 

To speed up the calculations, one can use the information on the Jacobian from 
previous steps and just iterate it. This is the idea of Broyden's method [3 1] and other 
gradient-free methods [32]. Given the Jacobian G„ for one step, the Jacobian for the 
next step G„+i is approximated by 

[5g„-G„(5x„] [g(x„+i) -x„+i] 
G„+i - G„ Ri J- — "2 O dx„ = — "2 dx„ (57) 

Here we have used the short notation (5x„ = x„_|_i— x„ and5g„ = g(x„+i)— g(x„), 
and (g) denotes the tensor product. The full Jacobian has to be calculated in the first 
step, and must possibly be updated after a number of steps. The method is a major 
improvement over Newton-Raphson, but it still requires the evaluation of at least one 
Jacobian. 

A very simple scheme which uses little computer time for one iteration is simple 
mixing. The [n + l)th guess x„+i is simply given by 

x„+i = x„ + A(g(x„) - x„) (58) 

with some appropriate A. Single iterations are thus cheap, but the method does not 
converge very well and many iteration steps are needed. The optimal values of A 
have to be established empirically (typically of the order of 0.1 or less). 

The simple mixing method can be improved by using information from previous 
steps to adjust A. For example, the results from the two previous steps can be used to 
determine a new value for A„ at every step [33] 
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\ — (Xn X„_i)^ ^^^^ 

y (g(x„) - x„ - g(x„_i) + x„_i)2 ■ 

A more systematic approach is "Anderson mixing", which is another standard 
method to solve systems of nonlinear equations in high dimensions [34, 35, 36J. 
From the remaining deviations d„ = g(x„) — x„ after n steps, one first determines 

Uij = (d„ - d„_i)(d„ - dn-j) (60) 
Vj = (d„ - d„_,)d„ (61) 

and then calculates 

x;^ = x„ + ^ Uij - 1 Vj (x„_i - x„) (62) 

i,i 

Sn = Si^n) + ^ C/.r'^J- (g(Xn-^) " g(x„)). (63) 

In the original Anderson mixing scheme, the {n + l)th guess x„-|_i is given by 
Xn+i = gn- More generally, the method can be combined with simple mixing 

x„+i = x;^ + A(g;^ - x;^). (64) 

This summarizes some methods which can be used to solve the self-consistent 
equations. Regardless of the particular choice, the single most time-consuming part 
in every SCF or EP calculation is the repeated solution of the diffusion equation 
(Eqs. (19) and (20)). Therefore the next task is to find an efficient method to solve 
equations of the type 

= ^Aq{r, t) - Wir)q{r, t). (65) 

with the initial condition q{r,0) = 1. Unfortunately, the two terms on the right hand 
side of this equation are incompatible. The Laplace operator A is best treated in 
Fourier space. The contribution of the field W{r) can be handled much more easily 
in real space. 

In SCF calculations of equilibrium mean-field structures, it is often convenient 
to operate in Fourier space, because crystal symmetries can be exploited. Relatively 
few Fourier components are often sufficient to characterize the structure of a phase 
satisfactorily. In Fourier space, Eq. (65) turns into a matrix equation 

^^ = -EAk,k'(iMk',i) (66) 
k,k' 

with 

Ak,k' = ^A:2,5k,k' + W{k, k') (67) 
and the initial condition q(k, 0) = (5(k, 0). The formal solution of this equation is 
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q{k,t) = ^[exp(-At)]k,k' g(k',0). 



(68) 



k' 



It can be evaluated by diagonalizing the matrix A. The approach has the advantage 
that the diffusion equation is solved exactly in time, /. e., one has no errors due to 
a time discretization. Unfortunately, the efficiency drops rapidly if the number of 
Fourier components and hence the dimension of the matrix A is very large. This 
hampers studies of fluctuations and dynamic phenomena, because one can no longer 
exploit crystal symmetries of the bulk structure, to reduce the number of variables. It 
is also an issue in SCF calculations if the synometry of the resulting structure is not 
yet known [37]. 

Therefore real space methods have received increased interest in recent years [37] . 
In real space, the number of degrees of freedom m can be increased more easily (the 
computing time typically scales like m or m In m). However, the time integration of 
the diffusion equation can no longer be performed exactly. 

Several different integration methods are available in the literature. Many of 
them approximate the effect of the Laplace operator by finite differences [31, 38]. 
Among the more elaborate schemes of this kind we mention the Crank-Nicholson 
scheme [31] and the Dufort-Frankel scheme [38]. 

In our own applications, we found that a different approach was far superior to 
finite difference methods: The pseudo- spectral split-operator scheme [39, 40]. This 
method combines Fourier space and real space steps, treating the contribution of the 
fields W{r) in real space, and the contribution of the Laplace operator in Fourier 
space. In that approach, the dynamical evolution during a single discrete time step of 
length h is formally described by 



where the effect of the exponential exp{hA) is evaluated in Fourier space. Taking as 
starting point a given q{r,t), the different steps of the algorithm are 

1. Calculate g'(r, t) = exp(— VK(r)/i/2) q{y, t) in real space for each r. 

2. Fourier transform (/'(r, i) (/'(k, 

3. Calculate g"(k,t) = exp(-/iZ\) g'(k,t) = exp(-/i(27rk/L)<^) 9'(k,i), where 
d is the dimension of space and L the system size. 

4. Inverse Fourier transform Q'"(k,t) — » q"{r,t) 

5. Calculate q{r, t + h) = exp{—W{r)h/2) q"{r, t) in real space for each r. 
5. Go back to 1. 

The algorithm thus constantly switches back and forth between the real space and 
the Fourier space. Even if this is done with efficient Fast Fourier Transform (FFT) 
routines, a single time step takes longer in this scheme than a single time step in one 
of the most efficient finite difference schemes. Moreover, the computing time now 
scales with m log m with the system size, instead of m. Nevertheless, we found that 
this is more than compensated by the vastly higher accuracy of the pseudo-spectral 
method, compared to finite difference methods. In order to calculate ^(r, t) with a 
given numerical accuracy, the time step h could be chosen an order of magnitude 
larger in the pseudo-spectral method than, e. g.,m the Dufort-Frankel scheme [41]. 



q{v, t + h)= e-^W''/^ e''^ g-^W'/^ g(r, t), 



(69) 
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3.5 

Extensions 

In order to make connection to experiments a variety of extension to the Gaussian 
chain model have been explored. Often it is desirable to include some of the 
architectural details for a better description on the length scale of the statistical 
segment length. Of course, the Gaussian chain model is only applicable if the 
conformational statistics on the length scale on which the composition varies is 
Gaussian. At strong segregation, ^ 1, or at surfaces to substrates or to 
vapor/vacuum, this assumption is often not fulfilled. 

The crossover from the Gaussian chain behavior on large length scales to the rod- 
like behavior on the length scale of the statistical segment length can be described by 
the worm-like chain model [ 1 7 ] . In this case, the end-segment distribution q depends 
both on the spatial coordinate r as well as on the orientation u, defined by the tangent 
vector of the space curve at position r [42, 43]. 

While the worm-like chain model still results in an analytical diffusion equation 
for the propagator, the structure of a polymer in experiments or computer simulation 
model often is not well-described by rod-like behavior. In order to incorporate 
arbitrary chain architecture and incorporate details of the molecular structure on 
all length scales, the sum (/ [R-]'Pi [R]) over the single chain conformations can 
be approximated by a partial enumeration over a large number of explicit chain 
conformations [18, 20, 21, 44, 45 J. Ideally they have been extracted from a computer 
simulation of the liquid phase, or they are obtained from a molecular modeling 
approach. Typically 10^ — 10"^ single chain conformations are employed in one- 
dimensional calculations. The longer the chains and the stronger the changes of the 
chain conformations in the spatially inhomogeneous field, the more conformations 
have to be considered. The sample size should be large enough that sufficiently 
many chains contribute significantly to the Boltzmann weight. The enumeration over 
the chain conformations is conveniently performed in parallel. To this end a small 
fraction of single chain conformations is assigned to each processor. Then, each 
processor calculates the Boltzmann weight of its conformations, the corresponding 
density profile, and the weight of its fraction of chain conformations. Subsequently, 
the total density profile is constructed by summing the weighted results of all 
processors. Typically, 64 or 128 processors are employed in parallel and a SCF 
calculation of a profile takes a few minutes on a CRAY T3E. As it is apparent, the 
detailed description of chain architecture comes at the expense of a large increase in 
computational demand. 

Also, in addition to the zero-ranged repulsion between unlike species, other 
interactions can be included into the theory. In order to capture some fluid-hke 
correlations and describe the details of packing in the vicinity of surfaces weighted 
density functionals have been used successfully. 

Electrostatic interactions have been included into the theory. This is important 
for tuning the orientation of self-assembled structures in block copolymers, and also 
to describe biological systems. 
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4 

Fluctuations 

4.1 

Examples of fluctuation effects and the role of JV' 

From the general expression of the free energy functional (11), we infer that the 
parameter combination fcsTV^ sets the scale of free energy fluctuations in a 
volume comparable to the size of the chain extension. Generally, if M is large, 
thermal fluctuations will not be important, and the self-consistent field theory will 
provide an adequate description. Qualitatively, the quantity M describes how strongly 
the molecules interdigitate: Fluctuation effects decrease with increasing number of 
interaction partners per molecule. There are, however, important exceptions where 
fluctuations change the qualitative behavior. Some examples of fluctuation effects in 
dense multicomponent noixtures shall be briefly mentioned: 

1) Unmixing transition in binary homopolymer blends 

In the vicinity of the critical point of a binary mixture one observes universal 
behavior, which mirrors the divergence of the correlation length of composition 
fluctuations. The universal behavior does not depend on the details of the 
system, but only on the dimensionality of space and the type of order parameter. 
Therefore binary polymer blends fall into the same universality class as mixtures 
of small molecules, metallic alloys, or the three-dimensional Ising model. In the 
vicinity of the critical point, Xc^ = 2 for a symmetric blend [14], the difference 
of the composition of the two coexisting phases - the order parameter m - 
vanishes like m ~ [xN — XcN)^, where the critical exponent /3 = 0.324 is 
characteristic of the 3D Ising universaUty class. This behavior is, however, only 
observable in the close vicinity of the critical point [46] 



of the mean-field theory, ^' < VGi ~ ^ [47, 44]. Outside 



the region marked by the Ginzburg criterion (70), fluctuations do not change 
the qualitative behavior. In the Umit of large AA, the region where fluctuations 
dominate the behavior becomes small. The cross-over from Ising to mean-field 
behavior as a function of chain length has attracted much experimental [48] and 
simulational interest [49, 50]. 

2) Fluctuations at the onset of microphase separation 

In case of diblock copolymers, the ordering is associated with composition 
waves of a finite periodicity in space. In contrast to homopolymer mixtures, 
not long-wavelength fluctuations are important, but those with wavevectors that 
correspond to the periodicity of the ordering. For symmetric block copolymers 



^ \XN - XcN] 
XcN 




(70) 



Fluctuations also shift the location of the critical point away from the prediction 
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one finds a fluctuation-induced first order transition (Brazovskii-mechanism [5 1]) 
and the shift of the transition temperature is [52] 

Xt^ = XtMF + 41.022Ar-V3 with xtMF = 10.495 (71) 



3) Capillary waves at the interface between two homopolymers 

If the systems contains interfaces, the local position of the interface can fluctuate. 
These capillary waves increase the area of the interface, and the free energy costs 
are proportional to the interface tension 7 [53]. Long-wavelength fluctuations 
are easily excited and influence measurements of the interfacial width w in 
experiments or computer simulations [54]. The resulting "apparent width" w 

depends on the lateral length scale L on which the interface is observed. Bq 
denotes a short-length scale cut-off on the order of the intrinsic width wq 
predicted by the self-consistent field theory [55] . Since is proportional to 
VM, the capillary wave contribution to the apparent interfacial width, measured 
in units of Rg, scales like 1 / vW [56]. 

4) Formation of microemulsion-like structure in the vicinity of Lifshitz points 

Capillary waves do not only broaden the width of the interface, but they can also 
destroy the orientational order in highly swollen lamellar phases. Those phases 
occur in mixtures of diblock-copolymers and homopolymers. The addition of 
homopolymers swells the distance between the lamellae, and the self-consistent 
field theory predicts that this distance diverges at Lifshitz points. However, 
general considerations show that mean-field approximations are bound to break 
down in the vicinity of Lifshitz points [58]. (The upper critical dimension is 
du = 8). This can be quantified by a Ginzburg criterion. Fluctuations are 
important if 

XtN 

The Ginzburg criterion has a simple heuristic interpretation: Fluctuations of 
the interface position destroy the long-ranged lamellar order. Those interface 
fluctuations are not Umited by the increase of the surface area (7 = 0), but rather 
by the bending stiffness k. The correlation length ^„ of the direction of a single, 
stiff interface is given by ^„ ~ exp{2TTK/kBT) [59]. If the stiffness k is of 
order unity, there is no long-range order of the normal direction of the interface, 
and the lamellar order will be replaced by a microemulsion-like structure. In 
the vicinity of a Lifshitz point self-consistent field theory predicts the interface 
stiffness to vanish like k kBT\fR which is compatible with the Ginzburg 
criterion (73) [60]. 
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Abbildung 1. Phase diagram for a ternary symmetrical A+B+AB blend from Monte Carlo 
simulations of the bond fluctuation model as a function of the incompatibility parameter xN 
and the homopolymer volume fraction $c- The chain length of all three species is equal, 
N = 32 corresponding to TV" = 91. 20 denotes a region of two-phase coexistence between an 
A-rich and a B-rich phase, 3 one of three-phase coexistence between an A-rich, a B-rich, and 
a lamellar phase. Slice through the three-dimensional system of linear dimension L = 4.7-Re 
and three periodic images are shown on the left hand side. From Ref. [57]. 



Similar arguments hold for a tricritical Lifshitz point, where a Lifshitz point 
happens to merge with a regular tricritical point. Here, fluctuations are important 
if 

^^xN-XtN^^^^y, (74) 
XtN 

Again this behavior is compatible with the behavior of the bending stiffness, 
K ~ ksTV^t^^^, as obtained fr om self-consistent field calculations [60] . 



5) Phase diagram of random copolymers 

IVIacrophase separation can also occur in random copolymers [61], which consist 
of a random sequence of Q blocks, each of which comprises either M segments 
of type A or M segments of type B. Macrophase separation occurs when xM is 
of order unity, /. e., independent from the number of blocks, and the A-rich and 
B-rich phases differ in their composition only by an amount of order 1 / y/Q. The 
strength of fluctuation effects can be quantified by the Ginzburg number [62] 

Gi ^ ^ with M - {p,nRli/Mf (75) 

where p,n = Ng denotes the monomer number density. While Gi decreases 
in the case of homopolymer mixtures like 1/JV, it increases quadratically in 
the number of blocks for random copolymers. In other words, if we increase 
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the number of blocks per molecule, fluctuation effects become stronger and 
stronger. Indeed, Monte Carlo simulations indicate the presence of a disordered, 
microemulsion-like structure in much of the parameter region where mean-field 
theory predicts macrophase separation. 

In all but the last example, fluctuation effects in three dimensions are controlled by 
the chain length N, in the sense that if A'' ~ AT — > oo, fluctuations effects can 
be made arbitrary small. This provides a formal limit in which self-consistent field 
theory becomes accurate. We emphasize, however, that in many practical examples 
fluctuations are very important. For instance, the broadening of the interface width 
by capillary waves is typically a sizeable fraction of the intrinsic width calculated by 
self-consistent field theory. 

In particle-based models, the mean-field approximation also neglects short- 
ranged correlations in the polymeric fluid, e. g., the fluid-like packing of the particles 
or the correlations due to the self-avoidance of the polymers on short length scales. 
There is no smaU parameter which controls the magnitude of these correlation 
effects, but they are incorporated into the coarse-grained parameters Re and x^- 



Gaussian fluctuations and Random Pliase Approximation 

As long as fluctuations are weak, they can be treated within the Gaussian approximation. 
One famous example of such a treatment is the random phase approximation (RPA), 
which describes Gaussian fluctuations in homogeneous phases. The RPA has been 
extended to inhomogeneous saddle points by Shi, Noolandi, and coworkers [63, 64]. 
In this section, we shall re-derive this generaUzed RPA theory within the formalism 
developed in the previous sections. 

Our starting point is the Hamiltonian ^J[f/, PF] of Eq. (11). We begin with 
expanding it up to second order about the saddle point G = G[U*,W*] of the SCF 
theory. To this end, we define the averaged single chain correlation functions [13, 63, 



with a,P = A OT B, where 0* is defined as in Eq. (28), and all derivatives are 
evaluated at the saddle point values of the fields, = {ill* + W*)/2 and = 
{ill* — W*)/2). As usual the sum J runs over the different types of polymers in 
the mixture, and nj denotes the total number of polymers of type J. In the case of a 
homopolymer mixture, the mixed correlation functions Kab{''^, r')> Kba{tc, r') are 
zero. 

With these definitions, the quadratic expansion of G[U, W] (11) can be written as 



4.2 



64] 




(76) 




(77) 
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with 5Wa = (iSU + SW)/2, SWb = {iSU - SW)/2. 

Since Q is an extremum, the Unear terms in SU and SW vanish. Next we expand 
this expression in SU and SW. To simplify the expressions, we follow Laradji et al., 
define [13, 64] 

E = Kaa + Kab + Kb A + Kbb (78) 
C = Kaa- Kab - Kb a + Kbb (79) 
A = Kaa + Kab - Kba - Kbb (80) 
A+ =Kaa-Kab + Kba-Kbb, (81) 
and adopt the matrix notation 



SA-K-SB = j d'^Y dVi^(r, r') SA{v) SB{r'). 



with the unity operator 1 = S{y—y'). Note that the operators S and C are symmetric, 
and A and A'^ are transposed to each other. In the new notation, Eq. (77) reads 

^ ' ^ SW-\-SW (82) 



VkBT{V/Ri) V 4xiV 

+ ^(SU-E-SU-SW-C-SW + iSW ■ A ■ SU + iSU ■ A+ ■ SW) . 

oV 

and can, after some algebra, be cast in the form 

= (83) 

y/MkBT{V/Ri) 

l^i^{SU-SU*)-S-{SU-SUn + ^SW-[l-^C]- SW 

with C = C-A-S-^-A+ and SU* = -iE''^ ■ A+ ■ SW. (84) 

The field U* + SU* corresponds to the saddle point of g[U, W] at fixed W = W* + 
SW. 

Both the fluctuations (SW^) and {{SU - SU*)) decrease Uke l/VA^. However, 
the fluctuations in both fields are qualitatively different. The operator S is positive, 
i. e., the eigenvalues of S are positive, therefore fluctuations in U remain small and 
the Gaussian approximation is justified in the limit of large V^. The operator 1 — 
XN/2C, on the other hand, may have negative eigenvalues at large x- In that case the 
saddle point under consideration is unstable and the true solution of the SCF theory 
corresponds to a different structure (phase). A saddle point which turns from being 
stable to unstable defines a spinodal or a continuous phase transition. In addition, 
the operator 1 — x-^/2C' also has zero eigenvalues in the presence of continuous 
symmetries. In all of these cases, the expansion for smaU SW becomes inaccurate 
and fluctuations in SW may give rise to qualitative deviations from the mean-field 
theory. 

Within the Gaussian approximation, the average {W{r)W{Y')) can readily be 
calculated: 
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{W{r)W{r')) = W*{r)W*{r') + ^^R] 



1 - 



(85) 



Using the relation between fluctuations of W and the composition, Eq. (26), we 
obtain 



{[Mr) 



6s(r)][^^(r') - Mr')]) - (Mr) - 4>BmMr') - Mr')) 
2Ri \^^_X^^ 



AfxN 
2Ri 



MxN 

1 



oo 



2 



c 



1 - 



XN^ 



XN^ 



(86) 



This is the general expression for Gaussian composition fluctuations in incompressible 
polymer blends derived from EP theory. The original derivation of Shi, Noolandi and 
coworkers [63, 64] uses as a starting point the density functional (53) and gives the 
identical result. A generalization for compressible blends can be found in Ref. [13]. 

The special situation of a homogeneous saddle point, corresponding to a homogeneous 
disordered phase, is particularly interesting. In that case, explicit analytical relations 
between the single chain partition function and the fields can be obtained, and 
one recovers the well-known random phase approximation (RPA). To illustrate this 
approach, we shall now derive the RPA structure factor for the case of a symmetrical 
binary homopolymer blend. 

Since the reference saddle point is spatially homogeneous, all correlation functions 
Kaf3{r, r') only depend on (r — r') and it is convenient to perform the calculations 
in Fourier space. We use the foUowing conventions for the Fourier expansion: 



(r) 



re 



(87) 



In homopolymer melts, the mixed single chain correlation functions Kab, Kba 

vanish, thus one has C = E = Kaa + Kbb and A = Z\+ = Kaa ~ Kbb- An 
explicit expression for the functional (/'^(r) = (j) A{(t> a{y)) in the limit of small fields 
Wa can be derived: 



<i>*A{v)=4>A\l-Y.M<\)WA^e^^' 



q#0 



= (PA 



gA{ci)WA{r')e 



iq(r— r') 



(88) 



A similar expressions holds for 4'*B{'r). Here g{q) = {2/x^){ex.p[—x] — 1 + x) with 
X = Rlq^ /6 denotes the Debye function. Hence, the Fourier transform of Kaa can 
be identified with the single chain structure factor: 
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Kaa{<^ = ^A5A(q), (89) 
and a similar expression holds Kbb- Using the definition of (7, Eq. (84), one obtains 



1 1 
+ 



Kaa K 



BB 



1 1 
+ 



(j^AQA <f>BgB 



Inserting this expression into Eq. (86), we can also obtain an explicit expression 
for the EP Hamiltonian n[W] = g[U* ,W] of binary blends within RPA: 



JfkBT{V/Ri) 8 ^0 4>a9a{c\) + ^s.gs(q) / 



(90) 



Combined with Eq. (26), this yields the well-known RPA expression for the 
collective structure factor in binary polymer blends 

4 1 1 

+ 7 TT - 2X^- (91) 



^/M{V|Ri){\4>A^-(t)B^?) 4>AgA{oL) ^b9b{ci) 

This calculation also justifies the use of Eq. (26) to calculate the fluctuations in the 
EP theory [28]. If we used the literal fluctuations in the EP theory according of 
Eq. (44), we would not recover the RPA expression but rather [28] 

(|<^Aq - <^BqP)EP 

- {\4>A^ ) y d r d r e ^ + 



J\f{V/Ri){l-2xNMBg{ci)) VJ\f{V/Ri) 

i\X X i2\ , (1 - #A0B)g(q) 
= (|^.,-0.,|)+ ^^^^^^^^ . (92) 

Both contributions are of the same order. 

This example illustrates how the RPA can be used to derive explicit, analytical 
expressions for Hamiltonians structure factors. The generalization to blends that also 
contain copolymers is straightforward. 



4.3 

Relation to Ginzburg-Landau models 

In order to make the connection to a Landau-Ginzburg theory for binary blends, we 
study the behavior of the structure factor at smaU wavevectors q for a synmietric 
mixture. Using g{q) Ri 1 + ''''^^ H we obtain: 

1 + 1 + ,93, 



VM{V/Ri){\(l}A^-(j>B^\') (I'A (t>B 18(/)A 18(As 
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If we assume composition fluctuations to be Gaussian, we can write down a free 
energy functional compatible with equation (91). 



'MkeTiV/Ri) v9v'(l//i?^)2(|0Aq-<^BqP) 
1/1 1 



2 V<^4 1 



2xN ) \cl>Ao - <^a|' (94) 

J2 d2 _2 d2 



2^V'/'A 1-(Aa 18(/)a 18(1 -(/-a) 

Note that this free energy functional is Gaussian in the Fourier coefficients of the 
composition, and, hence, the critical behavior still is of mean-field type. Transforming 
back from Fourier expansion for the spatial dependence to real space, we obtain for 
the free energy functional: 



'NkBT{VIRi) 




2 



360^(1 - (PA) 



2 dcl)\ V 27 360a(1-<^a) J 

(95) 



Expanding the free energy of the homogeneous system (Flory-Huggins free energy 
of mixing), we can restore higher order terms and obtain the Landau-de Gennes free 
energy functional for a symmetric binary polymer blend: 



360^(1 - (pA) J 

36(pA(l-(PA) ) 

A field-theory based on this simple expansion will already yield the correct three- 
dimensional Ising critical behavior. Note that the non-trivial critical behavior is 
related to higher order terms in W or (t>A — (f>B- Fluctuations in the incompressibUity 
field U or the total density 0a + </'b are not important. 



VAfkBT{V/Ri) 
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4.4 



Field-theoretic polymer simulations 

Studying fluctuations beyond the Gaussian approximation is difficult. Special types 
of fluctuations, e. g., capillary wave fluctuations of interfaces, can sometimes be 
described analyticaUy within suitable approximations [55]. The only truly general 
methods are however computer simulations. Here we shaU discuss two different 
approaches to simulating field theories for polymers: Langevin simulations and 
Monte Carlo simulations. 

4.4.1 

Langevin simulations 

As long as one is mainly interested in composition fluctuations (EP approximation, 

see section 3), the problem can be treated by simulation of a real Langevin process. 
The correct ensemble is reproduced by the dynamical equation 



where M(r, r') is an (arbitrary) kinetic coefficient, and 9{r, t) is a stochastic noise. 
The first two moments of 6 are fixed by the fluctuation-dissipation theorem 



The choice of M(r,r') determines the dynamical properties of the system. For 
example, Af (r, r') = 5{y — r')/TkBT corresponds to a non-conserved field, while 
field conservation can be enforced by using a kinetic coefficient of the form M (r, r') = 
Vr^(r — r')Vr'. Different forms for the Onsager coefficient A will be discussed 
in Sec. 5. In each step of the Langevin simulation one updates all field variables 
simultaneously [28] and the the self-consistent equations for the saddle point U* (W) 
have to be solved. 

The approach is commonly referred to as external potential dynamics (EPD). 
A related approach has originally been introduced by Mauri ts and Fraaije [30]. 
However, these authors do not determine U*{W) exactly, but only approximately 
by solving separate Langevin equations for real fields Wa and Wb- This amounts to 
introducing a separate Langevin equation for a real field iU {i. e., an imaginary U in 
our notation) in addition to Eq. (97). 

As the fluctuation effects discussed earlier (cf. Sec. 4.1) are essentially caused 
by composition fluctuations, the EP approximation seems reasonable. Nevertheless, 
one would also like to study the full fluctuating field theory by computer simulations. 
In attempting to do so, however, one faces a serious problem: Since the field iU is 
imaginary, the Hamiltonian Q is in general complex, and the "weight" factor cxp[— fj] 
can no longer be used to generate a probability function. This is an example of a 
sign problem, as is weU-known from other areas of physics [65], e. g., lattice gauge 
theories and correlated fermion theories. A number of methods have been devised to 




(97) 



W=0 



{0{r, t)e{r', t)) = 2kBTM{r, r') d{r - r') 6{t - f'). 



(98) 
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handle complex Hamiltonians or complex actions [66, 67, 68, 69, 70]. Unfortunately, 
none of them is as universally powerful as the methods that sample real actions 
(Langevin simulations, Monte Carlo, etc. ). 

Fredrickson and coworkers [71, 72, 73, 74] have recently introduced the complex 
Langevin method [66, 67] into the field of polymer science. The idea of this method 
is simply to extend the real Langevin formalism, e. g., as used in EPD, to the case 
of a complex action. The drift term dH/5W in Eq. 97 is replaced by complex drift 
terms SQ[U,W]/SU and SG[U,W]/5W, which govern the dynamics of complex 
fields U and W. This generates a diffusion process in the entire complex plane for 
both U and W. One might wonder why such a process should sample line integrals 
over U and W. To understand that, one must recall that the integration paths of 
the line integrals can be distorted arbitrarily in the complex plane, as long as no 
pole is crossed, without changing the result. Hence a complex Langevin trajectory 
samples an ensemble of possible line integrals. Under certain conditions, the density 
distribution converges towards a stationary distribution which indeed reproduces the 
desired expectation values [75]. Unfortunately, these conditions are not generally 
satisfied, and cases have been reported where a complex Langevin process did not 
converge at all, or converged towards the wrong distribution [68, 76]. Thus the 
theoretical foundations of this methods are still under debate. In the context of 
polymer simulations, however, no problems with the convergence or the uniqueness 
of the solution were reported. 

In our case, the complex Langevin equations that simulate the Hamiltonian (11) 
read 

dU{r,t) Sg[U,W]/kBT 



dt SU{r,t) 
dW{r,t) 6g[U,W]/kBT 



Ou{v,t) (99) 
+ 6w{v,t), (100) 



dt 5W{v,t) 

where 6u (r, t), 6w (r, t) are real Gaussian white noises which satisfy the fluctuation 
dissipation theorem 



(61) =0 {e{Y,t)e{T',t))=2 5{T-T')5{t-t'). (101) 
The partial derivatives are given by 

(102) 



5g[U,W\/kBT _ 
SW{r,t) ~ 2Ri 

Sg[U,W]/kBT _ 



w 



SU{r,t) 2i?f 



1 



(103) 



with and (f)*^ defined as in Eq. (28). 

If the noise term is turned off, the system is driven towards the nearest saddle 
point. Therefore, the same set of equations can be used to find and test mean- 
field solutions. The complex Langevin method was first apphed to dense melts of 
copolymers [71], and later to mixtures of homopolymers and copolymers [77] and to 
diluted polymers confined in a slit under good solvent conditions [74]. Fig. 2 shows 
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(a) (b) (c) 

Abbildung 2. Averaged densities across tlie order-disorder transition in a two-dimensional 
ternary system with A,B homopolymers and A-B copolymers (20 % homopolymer volume 
fraction), as obtained from Complex Langevin simulation runs. 

examples of average "density" configurations (^*) for a ternary block copolymer/ 
homopolymer system above and below the order/disorder transition. 



4.4.2 

Field-theoretic Monte Cario 

As an alternative approach to sampling the fluctuating field theory, Diichs et al. [41, 
77] have proposed a Monte Carlo method. Since the weight exp[— ^] is not positive 
semidefinite, the Monte Carlo algorithm cannot be applied directly. To avoid this 
problem, we split Q into a real and an imaginary part C/^ and iQ^ and sample only 
the real contribution exp(— [68] . The imaginary contribution is incorporated into 
a complex reweighting factor cxp(— i^?^). Statistical averages over configurations j 
then have to be computed according to 

{A) = % (104) 

i. e., every configuration is weighted with this factor. Furthermore, we premise that 
the (real) saddle point iU* contributes substantially to the integral over the imaginary 
field i U and shift the integration path such that it passes through the saddle point. The 
field iU is then represented as 

iU = iU*[W]+iw, (105) 

where u) is real. 

The Monte Carlo simulation includes two different types of moves: trial moves 
of W and w. The moves in ui are straightforward. Implementing the moves in W is 
more involved. Every time that W is changed, the new saddle point iU* [W] must be 
evaluated. In an incompressible blend, the set of self consistent equations 

'^A(r)+.^s(r) = l, Vr (106) 



must be solved for iU*. Fortunately, this does not require too many iterations, if W 
does not vary very much from one Monte Carlo step to another. 



Fluctuations and dynamics in self-consistent field theories 



33 



Compared to the Complex Langevin method, the Monte Carlo method has 
the advantage of being well founded theoretically. However, it can become very 
inefficient when spreads over a wide range and the reweighting factor oscillates 
strongly. In practice, it relies on the fact that the integral is indeed dominated by one 
(or several) saddle points. 

Can we expect this to be the case here? To estimate the range of the reweighting 
factor, we briefly re-inspect the Hubbard-Stratonovich transformation of the total 
density (f>A+4>B that lead to the fluctuating field U. For simplicity we consider a one- 
component system. In a compressible polymer solution or blend, the contribution of 
the repulsive interaction energy to the partition function can be written as 



exp 



qN 



(107) 



(with ^0 = 1 in a melt, and = in a good solvent). The Hubbard-Stratonovich 
transformation of this expression is proportional to 



J VUe^p[-g{^ Jd^'vU^ + iJ d^r ^{4>a + 4>b - M} 



(108) 



Thus U should to be distributed around the saddle point with a width proportional 
to K. The method should work best for very compressible solutions with small k. 
In contrast, in a compressible blend, the contribution (107) is replaced by a delta 
function constraint 6{(j)A + (I>b — 1). The fluctuating field representation of this 
constraint is 

j VUeM-iQ j d'^r + ^b - 1)]- (109) 

The only forces driving U towards the saddle point are now those related to the 
single chain fluctuations (cf. Sec. 4.2), and U will be widely distributed on the 
imaginary axis. Thus it is not clear whether the Monte Carlo method will work for 
incompressible systems. 

Indeed, simulations of an incompressible ternary A+B+AB homopolymer/ copolymer 
blend showed that the values of the argument of the reweighting factor cover a 
wide range (Fig. 3). Hence computing statistical averages becomes very difficult, 
because both the numerator and the denominator in Eq. (104) are subject to very 
large relative errors. 

Fortunately, a closer inspection of the data reveals that the reweighting factor 
is entirely uncorrelated with the fluctuating field W. Since the latter determines all 
quantities of interest to us (Eq. (25) and (26)), the averages can be computed without 
reweighting. Moreover, the time scales for variations of W and uj are decoupled. By 
samphng u) more often than W, they can be chosen such that the time scale of ui is 
much shorter than that of W. Finally, the dynamics of the actual simulation, which 
is governed by the real factor cxp(—G^), is virtually the same as that of a reference 
simulation with lu switched off [77]. 

Combining all these observations, we infer that the fluctuations in w do not 
influence the composition structure of the blend substantially. A similar decoupling 
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1000 2000 3000 4000 5000 

MCS 

Abbildung 3. Argument of the reweighting factor in a simulation of an incompressible 
blend. From Ref. [77]. 




Abbildung 4. Snapshots of (j)*^ for the same ternary A+B+AB system as in Fig. 2, at 
70 % homopolymer volume fraction, for different above and below the order-disorder 
transition, as obtained from Monte Carlo simulations with lj switched off (EP theory). From 
Ref. [80]. 

between composition and density fluctuations is suggested by other studies [78, 
79, 44]. The composition structure can equivalently be studied in a Monte Carlo 
simulation which samples only W and sets u; = 0. We have demonstrated this for 
one specific case of an incompressible blend and suspect that it may be a feature 
of incompressible blends in general. The observation that the fluctuations in W 
and Lu are not correlated with each other is presumably related to the fact that the 
(vanishing) density fluctuations do not influence the composition fluctuations. If 
that is true, we can conclude that field-theoretic Monte Carlo can be used to study 
fluctuations in polymer mixtures in the limits of high and low compressibilities. 
Whether it can also be applied at intermediate compressibilities will have to be 
explored in the future. 

Fig. 4 shows examples of snapshots for the ternary system in the vicinity of the 
order-disorder transition. 

We should note that the Monte-Carlo simulation with u; = effectively samples 
the EP Hamiltonian. This version of field-theoretic Monte Carlo is equivalent to 
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the real Langevin method (EPD), and can be used as an alternative. Monte Carlo 
methods are more versatile than Langevin methods, because an almost unUmited 
number of moves can be invented and implemented. In our applications, the W and 
w-moves simply consisted of random increments of the local field values, within 
ranges that were chosen such that the MetropoUs acceptance rate was roughly 35%. 
On principle, much more sophisticated moves are conceivable, e. g., collective moves 
or combined EPD/Monte Carlo moves (hybrid moves [81]). On the other hand, EPD 
is clearly superior to Monte Carlo when dynamical properties are studied. This will 
be discussed in the next section. 



Dynamics 

5.1 

Onsager coefficients and dynamic SCF theory (DSCFT) 

In order to describe the diffusive dynamics of composition fluctuations in binary 
mixtures one can extend the time-dependent Ginzburg Landau methods to the free 
energy functional of the SCF theory. The approach relies on two ingredients: a 
free energy functional that accurately describes the chemical potential of a spatially 
inhomogeneous composition distribution out of equiUbrium and an Onsager coefficient 
that relates the variation of the chemical potential to the current of the composition. 

We imphcitly assume that the out of equihbrium configuration can be described 
by the spatially varying composition, but that we do not have to specify the details 
of the molecular conformations. This can be justified if the time scale on which the 
composition changes is much larger than the molecular relaxation time, such that 
the chain conformations are always in "local equilibrium" with the instantaneous 
composition. The time scale for the chain conformations to equilibrate is set by 
T = Rl/D, where D denotes the single chain self-diffusion constant. Moreover, we 
assume that the free energy obtained from the SCF theory also provides an accurate 
description even out of equilibrium, /. e., not close to a saddle point. 

Then, one can relate the spatial variation of the chemical potential to a current of 
the composition. 



J^{r,t) denotes the current density at position r at time t. There is only one 
independent chemical potential, because of the incompressibUity constraint and 
it is given by /x^ = with J"* [c^a] from Eq. (53). 



5 




(110) 



V^keT 



(111) 



RPA 



ixN - 2)(1 - 2Mr)) - -f^Mr) 



(112) 
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where if2^ and are real fields. In the last equation we have used the exphcit 
expression for the chemical potential within RPA (cf. Eq. (96)) for a symmetric 
binary polymer blend. 

The kinetic Onsager coefficient yl0(r, r') relates the "chemical force" acting 
on a monomer at position r' due to the gradient of the chemical potential to 
the concomitant current density at position r. This describes a purely relaxational 
dynamics, inertial effects are not captured. can be modeled in different ways: The 
simplest approach would be local coupling which results in the Onsager coefficient 



~ -52^^ ,t)(t>B[r ,t)Rl5[T - r ) (local) (113) 



Rl Rl 

The composition dependence accounts for the fact that currents of A- and B- 
densities have to exactly cancel in order to fulfill the incompressibility constraint. 
This local Onsager coefficient completely neglects the propagation of forces along 
the backbone of the chain and monomers move independently. Such a local Onsager 
coefficient is often used in calculations of dynamic models based on Ginzburg- 
Landau type energy functionals for reasons of simplicity [82, 83, 84]. 

Bearing in mind that the connectivity of monomeric units along the backbone of 
the polymer is an essential ingredient of the single chain dynamics it is clear that non- 
local coupUng should lead to a better description. In the Rouse model forces acting 
on a monomer caused by the chain cormectivity are also taken into account [85, 86]. 
This leads to a kinetic factor that is proportional to the intramolecular pair-correlation 
function [87, 88, 89], P(r, r') 



= ^-^^(r '*)<P-B(r >*)-P(r>r ) (Rouse). (114) 



Rl Rl 

The difficulty in using Rouse dynamics Ues in the computational expense of 
calculating the pair-correlation function in a spatially inhomogeneous environment 
at each time step. If the system is still fairly homogeneous (e. g., in the early stages 
of a demixing process), one can use the pair-correlation function of a homogeneous 
melt, as it is given through RPA [16]. This leads to the following Onsager coefficient 
in Fourier space: 



Anouseiq) = Aioc^i{q)g{q) = -^<i)A<pB- — ^ (115) 



X is defined as a; = q^Rl/6 and g is the Debye function. Another model for non- 
local coupling is the reptation model [86, 90] which is appropriate for polymer melts 
with very long chains, i. e., entangled chains. 

Since the total amount of A-component is conserved, the continuity equation 
relates the current of the composition to its time evolution: 

'-^+VMr,t) = (116) 
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Eqs. (116) and (110) lead to the following diffusion equation: 

d(j)A{r,t) 



dt 



Vr J dVyl^(r, r')Vr'M<^(r', t) + r]{r, t) 



(117) 



the last term representing noise that obeys the fluctuation-dissipation theorem. After 
Fourier transformation this diffusion equation adopts a very simple form: 



d(j)A{(l, t) 

dt 



RPA VMkBTA^{qR^f 



RPA 



2{xN-2) + -{qR,f 



1 



+ 



Ag{q) <i)B9{q) 



-2xN 



(118) 
6A + r?(119) 

(120) 



The last two equations are the explicit expressions within RPA, and g{q) denotes the 
debye function. 

We have now found all necessary equations to numerically calculate the time 
evolution of the densities in a binary polymer mixture. This leads us the following 
procedure to which we refer as the dynamic self consistent field theory (DSCFT) 
method: 

(1) Calculate the real fields iQ\ and if2% that "create" the density distribution 
according to Eq. (52) using the Newton-Broyden scheme. 

(2) Calculate the chemical potential according to Eq. (111). 

(3) Propagate the density in time according to Eq. (118) using a simplified Runge- 
Kutta method. 

(4) Go back to (1). 



5.2 

External Potential Dynamics (EPD) 

Instead of propagating the composition in time, we can study the time evolution of 
the exchange potential W. In equilibrium the density variable (I)a — 4'b and the field 
variable W = iQa — i^B are related to each other via cjyA — 4'b = —W/xN, see 
Eq. (25). We also use this identification to relate the time evolution of the field W 
to the time dependence of the composition. Since the composition is a conserved 
quantity and it is linearly related to the field variable, we also expect the field W 
with which we are now describing our system to be conserved. Therefore we can 
use the free energy functional T-L\W] from Eq. (42) and describe the dynamics of 
the field W through the relaxational dynamics of a model B system, referring to the 
classification introduced by Hohenberg and Halperin [91]. 



dW{v,t) 
di 



Vr j d'^r ylEPD(r, r')Vr'm(r', t) + r?(r, t) (121) 
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with the chemical potential fiyy being the first derivative of the free energy Ti[W] 
with respect to the order parameter W, 



l^w jr) ^ 6H[W{r)] ^ VMksTW + xN [^^(r) - (/>^(r)] 
keT 5W{r) Ri 2xN ' ^ ' 

^EPD is a kinetic Onsager coefficient, and r] denotes noise that satisfies the fluctuation- 
dissipation theorem. The Fourier transform of this new diffusion equation is simply: 

dW{q,t) _ VMkBTAEPD{qRe)^W + xN{rA-rB)^^ (^23) 



dt 2xN 

RPA VMksTAEPDiqRef W [l - 2xN4>aM<i)) 



+ V (124) 



Ri+^ 2xN 

ri{q,t) is white noise that obeys the fluctuation-dissipation theorem. In the last 
equation we have used the Random-Phase Approximation for the Hamiltonian of the 
EP theory (cf. Eq. (90)). We refer to the method which uses this diffusion equation 
as the external potential dynamics (EPD) [30]. As shown in Sec. 4.4.1, any Onsager 
coefficient (in junction with the fluctuation-dissipation theorem) will reproduce the 
correct thermodynamic equilibrium. But how is the dynamics of the field W related 
to the collective dynamics of composition fluctuations? 

Comparing the diffusion equations of the dynamic SCF theory and the EP 
Dynamics, Eqs. (120) and (124), and using the relation W{q,t) = —2xN<j)A{q,t), 
we obtain a relation between the Onsager coefficients within RPA: 

vlEPD(g) = - ^ A^{q) (125) 

In particular, the non-local Onsager coefficient ylRouse that mimics Rouse-like dynamics 
(cf . Eq. (1 14) corresponds to a local Onsager coefficient in the EP Dynamics 

^ksTA^^^.^^^Uq) ^ ^^^^ (^26) 

Generally, one can approximately relate the time evolution of the field w to the 
dynamic SCF theory [30]. The saddle point approximation in the external fields 
leads to a bijective relation between the external fields if^% and (j)A, <t^B- In 
the vicinity of the saddle point, we can therefore choose with which of the two sets 
of variables we would like to describe the system. Using the approximation, 

VrP(r,r')^-Vr'P(r,r') (127) 

one can derive Eq. (126) without invoking the random phase approximation. The 
approximation (127) is obviously exactly valid for a homogeneous system, because 
the pair-correlation function only depends on the distance |r— r'| between two points. 
It is expected to be reasonably good for weakly perturbed chain conformations. 
To study the time evolution in the EPD we use the following scheme [28]: 
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(0) Find initial fields, W = if2^- ifi% and U = ifi'^ + i^B' ^^^^ create the initial 
density distribution 

(1) Calculate the chemical potential ij,w according to Eq. (122). 

(2) Propagate W according to Eq. (123) using a simplified Runge-Kutta method. 

(3) Adjust U to make sure the incompressibiUty constraint (/)^ + 0^ = 1 is fulfilled 
again using the Newton-Broyden method. 

(4) Go back to (1). 

The EPD method has two main advantages compared to DSCFT: First of all it 
incorporates non-local coupUng corresponding to the Rouse dynamics via a local 
Onsager coefficient. Secondly it proves to be computationally faster by up to one 
order of magnitude. There are two main reasons for this huge speed up: In EPD 
the number of equations that have to be solved via the Newton-Broyden method to 
fulfill incompressibility is just the number of Fourier functions used. The number of 
equations in DSCFT that have to be solved to find the new fields after integrating the 
densities is twice as large. On the other hand, comparing the diffusion equation (117) 
used in the DSCFT method with equation (121) in EPD, it is easily seen, that the 
right hand side of the latter is a simple multiplication with the squared wave vector 
of the relevant mode, whereas the right hand side of equation (117) is a complicated 
multiplication of three spatially dependent variables. 

6 

Extensions 

The dynamic self-consistent field theory has been widely used in form of MESODYN [92]. 
This scheme has been extended to study the effect of shear on phase separation or 
microstructure formation, and to investigate the morphologies of block copolymers 
in thin films. In many practical applications, however, rather severe numerical 
approximations (e. g., very large discretization in space or contour length) have been 
invoked, that make a quantitative comparison to the the original model of the SCF 
theory difficult, and only the qualitative behavior could be captured. 

More recently, stress and strain have been incorporated as conjugated pair of 
slow dynamic variables to extend the model of the SCF theory. This allows to 
capture some effects of viscoelasticity [73]. Similar to the evaluation of the single 
chain partition function by enumeration of explicit chain conformations, one can 
simulate an ensemble of mutually non-interacting chains exposed to the effective, 
self-consistent fields, U and W, in order to obtain the densities [93]. Possibilities to 
extend this scheme to incorporate non-equilibrium chain conformations have been 
explored [94]. 
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Abbildung 5. Mean-field phase diagram for a ternary symmetrical A+B+AB blend as a 
function of the incompatibility parameter xN and the homopolymer volume fraction cPh- The 
chain lengths of the copolymers is five time that of homopolymers. 2<j) denotes a region of two- 
phase coexistence between an A-rich and a B-rich phase, 3 4> one of three-phase coexistence 
between an A-rich, a B-rich, and a lamellar phase. From Ref. [77]. 

7 

Applications 

7.1 

Homopolymer-copolymer mixtures 

As discussed in Sec. 4.1, one prominent example of a situation where the SCF theory 
fails on a qualitative level is the microemulsion channel in ternary mixtures of A and 
B homopolymers and AB diblock copolymers. Fig. 5 shows an example of a mean- 
field phase diagram for such a system. Four different phases are found: A disordered 
phase, an ordered (lamellar) phase (see Fig. 4, right snapshot), an A-rich and a B- 
rich phase. The SCF theory predicts the existence of a point where all three phases 
meet and the distance of the lamellar sheets approaches infinity, an isotropic Lifshitz 
point [95, 96]. 

It seems plausible that fluctuations affect the Lifshitz point. If the lamellar 
distance is large enough that the interfaces between A and B sheets can bend around, 
the lamellae may rupture and form a globally disordered structure. A Ginzburg 
analysis reveals that the upper critical dimension of isotropic Lifshitz points is as high 
as 8 (see also Sec. 4.1). The lower critical dimension of isotropic Lifshitz points is 
d* = 4 [97, 98]. Thus fluctuations destroy the transition in two and three dimensions. 

Indeed, the experimentally observed phase behavior differs substantially from 
the mean-field phase diagram. An example is shown in Fig. 6. The Lifshitz point 
is destroyed, the three phase coexistence region between the lameUar phase and the 
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Abbildung 6. Experimental phase diagram of a PDMS-PE/PDMS/PE blend with 
copolymers about five times as long as the homopolymers. From Ref. [99]. Reproduced by 
permission of the Royal Society of Chemistry. 

demixed A-rich and B-rich phases is removed and gives way to a macroscopically 
mixed phase. A transmission electron micrograph of this phase is shown and compared 
with the structure of the lamellar phase in Fig. 7. 

In order to study this effect, Diichs et al. have performed field-theoretic Monte 
Carlo simulations of the system of Fig. 5 [77, 80, 100], in two dimensions. (For 
the reasons explained in Sec. 4.4.2, most of these simulations were carried out in 
the EP approximation). Some characteristic snapshots were already shown in Fig. 4. 
Here we show another series of snapshot at xTV = 12.5 for increasing homopolymer 
concentrations (Fig. 8). For all these points, the self-consistent field theory would 
predict an ordered lamellar phase. In the simulation, only configurations with lower 
homopolymer concentrations exhibit (defective) lamellar order. At higher homopolymer 
concentration, the order breaks up and a microemulsion emerges. 

The phase transition can be identified by defining appropriate order parameters. 
We use anisotropy parameters, which are defined in terms of the Fourier transform 
-F(q) of images as those shown in Fig. 8: 



Fn (q) is zero for isotropic (disordered) configurations and nonzero for anisotropic 
(lamellar) configurations. The anisotropy information can be condensed into a single 
dimensionless parameter 




(128) 
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K := (129) 



with the normalization 

(130) 



Examples of anisotropy parameters are shown in Fig. 9(left). They are clearly suited 
to characterize the phase transition between the disordered phase at low xiV and the 
lamellar phase at high x^^- For comparison, the same system was also examined with 
complex Langevin simulations. The results are shown in Fig. 9(right). The transition 
points are the same. 

Unfortunately, the order of the transition could not yet determined from these 
simulations. In the self-consistent field theory, the transition is continuous. According 
to theoretical predictions [51, 52], fluctuations shift it to lower xiV and change the 
order of the transition. Our simulation data are consistent both with the assumption 
that the transition is continuous or weakly first order In particular, we have established 
that the miscibihty gap between the ordered and the disordered phase is very small 
(data not shown). 

Besides the order/disorder transition, we also need to determine the demixing 
transition. This was done in the grand canonical ensemble by monitoring the difference 
between A and B monomer densities. It exhibits a jump at the phase transition. All 
phase transition points are summarized in the phase diagram of Fig. 10. It is in good 
qualitative agreement with the experimental phase diagram of Fig. 6. In particular, 
we reproduce the cusp-like region of the microemulsion. 




Abbildung 7. Transmission electron micrographs from symmetric PE/PEP/PE-PEP blends: 
(a) fluctuating lamellae (b) microemulsion phase. From Ref. [99]. Reproduced by permission 
of the Royal Society of Chemistry. 
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Abbildung 8. Snapshots of for the ternary A+B+AB system of Fig. 5 at xN = 12.5 for 
homopolymer volume fraction (j)H = 0.74,0.78,0.8, 0.82 (first row from left), and 0.84, 0.86, 
and 0.9 (second row from left). 




Abbildung 9. Anisotropy parameters F2 (solid) and F4, (dashed) vs. for different 
homopolymer volume fractions at the order/disorder transition. Left: From Monte Carlo 
Simulations (EP theory) Right: From Complex Langevin simulations. After Ref. [77] 

The simulations allow to analyze the phase transition and the microemulsion 
phase in more detail. Fig. 1 1 compares different characteristic length scales of the 
system at the order/disorder transition. One length scale, io, is obtained from the 
maximum of ^o(<z) (Eq. (128)) and gives the typical distance between layers. The 
other is the mean absolute curvature radius Dc, and characterizes the length scale for 
layer bending. Whereas the typical layer distance increases with the homopolymer 
concentration, the curvature radius remains roughly constant. The microemulsion 
phase transition takes place at the value of (f)H where the two lines cross. Thus 
we can identify the mechanism by which fluctuations generate the microemulsion: 
The lamellae break up when their width grows larger than the length scale which 
characteristizes the boundary fluctuations. 
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Abbildung 10. Phase diagram from simulations for the same system as above. The circles 
show the results from Monte Carlo simulations, the squares those from the complex Langevin 
simulations. The dotted lines correspond to the mean-field prediction from Fig. 5. From 
Ref. [77]. 



To summarize, the example of homopolymer/copolymer mixtures demonstrates 
nicely how field-theoretic simulations can be used to study non-trivial fluctuation 
effects in polymer blends within the Gaussian chain model. The main advantage of 
these simulations is that they can be combined in a natural way with standard self- 
consistent field calculations. As mentioned earlier, the self-consistent field theory is 
one of the most powerful method for the theoretical description of polymer blends, 
and it is often accurate on a quantitative level. In many regions of the parameter 
space, fluctuations are irrelevant for large chain lengths (large Af) and simulations 
are not necessary. Field-theoretic simulations are well suited to complement self- 
consistent field theories in those parameter regions where fluctuation effects become 
important. 



7.2 

Spinodal decomposition 

To illustrate the application of the method to study the dynamics of collective 
composition fluctuations, we now consider the spontaneous phase separation e., 
spinodal decomposition) that ensues after a quench into the miscibihty gap at the 
critical composition, 4' A = 4>b = 1/2, of the symmetric blend. Different time 
regimes can be distinguished: During the early stages, composition fluctuations 
are amplified. Fourier modes with a wavevector q smaller than a critical value qc 
exponentially grow and modes with different wavevectors are decoupled. As the 
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Abbildung 11. Characteristic lengtlr scales at ~ 12-5 as a function of liomopolymer 
concentration 4>h- Dc (circles) shows the mean (absolute) curvature radius, and Lq the 
preferential length scale of layer distances, obtained from the maximimi of Fo (g) (Eq. (128)). 
FromRef. [77] 



amplitude of the modes continues to grow, the composition in the domains begins to 
saturate and different modes begin to interact. At even later stages of phase ordering, 
hydrodynamics dominates the coarsening. 

In the following, we restrict our attention to the early stages of spinodal decomposition. 
In the analysis of experiments one often uses the Landau-de Gennes functional 
(96). The Cahn-Hilliard-Cook theory [101] for the early stages of phase separation. 
This treatment predicts that Fourier modes of the composition independently evolve 
and increase exponentially in time with a wavevector-dependent rate, (I>a{q, i)^ ~ 
cxp[R{q)t]. Therefore, it is beneficial to expand the spatial dependence of the 
composition in our dynamic SCF or EP calculations in a Fourier basis of plane waves. 
As the Unearized theory suggests a decoupling of the Fourier modes at early stages, 
we can describe our system by a rather small number of Fourier modes. 

Using Eq. (120) one reads off the rate R{q): 



i?(q) = 4 



UksTA^iqEl) 



r: 



■d+2 



xiV-2. 



and obtains for local dynamics and Rouse-like dynamics 



D 
D 



= {qR,f{xN-2) 
= {qR,f{xN-2)9{q) 



(131) 

(132) 
(133) 
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Abbildung 12. Growth rates of density modes for a quench of a symmetric mixture to 
XN = 5. Below the critical wave vector qc the density modes increase spontaneously. Modes 
with larger wave vectors are damped. From Ref. [28]. 



respectively, where the critical wavevector is given by qcR = \/%xN — 2). During 
the early stage, composition modes with wavevectors q < Qc will grow. In the case 
of local dynamics, the fastest growth occurs at qc/ \/2. In the case of Rouse-like 
dynamics, the maximum of the growth rate is shifted to smaller wavevectors. The 
expressions above use the RPA which becomes accurate at weak segregation (WSL). 
Of course, the applicability of the linearized theory of spinodal decomposition using 
the Landau-de Gennes free energy functional is severely limited: In the ultimate 
vicinity of the critical point, the free energy functional does not even give an 
accurate description of the thermodynamic equilibrium, and the relaxation times are 
characterized by a non-classical, dynamical critical exponent. 

In the opposite limit of strong segregation (SSL) 2, the free energy can 

also be expressed in terms of a square gradient functional, and one obtains qcR = 
\/6{x^ " 2). At intermediate segregation, which is most relevant to experiments, 
the square-gradient approximation breaks down and there is no straightforward way 
to interpolate between the simple analytical expression in the weak and strong 
segregation hmit. 

The predictions of the rate R{q) for the weak and the strong segregation Umit 
are presented in Fig. 12. As expected the results of the dynamic SCF theory with a 
local Onsager coefficient fall between the results of the Cahn-Hilliard-Cook theory 
in the two Umits of weak and strong segregation. Comparing the predictions for the 
kinetics of phase separation to experiments or simulations, it often remains unclear 
whether deviations are caused by approximating the free energy functional by the 
Landau-de Germes expressions or the Onsager coefficient. Indeed, earUer Monte 
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Abbildung 13. Global structure factor versus wave vector for different times for a quench 
from the mixed state at ~ 0.314 to ~ 5- Lines represent Monte Carlo results, 
filled symbols dynamic SCF theory results. Panel (a) compares dynamic SCF theory using 
a local Onsager coefficient with the Monte Carlo simulations. Local dynamics obviously 
overestimates the growth rate and shifts the wavevector that corresponds to maximal growth 
rate to larger values. Panel (b) compares dynamic SCF theory using a non-local Onsager 
coefficient that mimics Rouses dynamics with Monte Carlo results showing better agreement. 
FromRef. [28]. 



Carlo simulations found rather pronounced deviations [102] from the prediction of 
the Cahn-HilHard-Cook theory [101]. 

To overcome this difficulty we use EP dynamics or dynamic SCF theory, which 
do not invoke a gradient expansion of the free energy functional and yield an accurate 
description of the free energy costs of composition fluctuations. Quantitatively 
comparing dynamical calculations to computer simulations of the bond fluctuation 
model [ 1 03] , we investigate the relation between the dynamics of collective composition 
fluctuations and the underlaying dynamics of single chains, which is encoded in the 
Onsager coefficient A. The parameters of the bond fluctuation model can be related 
to the coarse-grained parameters, xiV, Re and Af, of the standard SCF model and the 
static properties of the bond fluctuation model (e. g., the bulk phase behavior [44, 50] 
and interface properties [43, 45, 54, 55]) have been quantitatively compared to SCF 
theory. Note that the mapping between the particle-based simulation model and the 
SCF model does not involve any adjustable parameter: the Flory-Huggins parameter 
XN is identified via the energy of mixing, the length scale is set by and the 
time scale is determined by the self-diffusion constant of a single polymer chain 
in a dense melt. All these quantities are readily measurable in the simulations. In 
the following example, we consider a quench to = 5- This is a compromise: 
On the one hand it is sufficiently below the critical point, XcN — 2, such that 
critical composition fluctuations are not important. On the other hand the interfaces 
between the coexisting phases are still wide enough such that the polymers are 
still describable by the Gaussian chain model on the length scale of the (intrinsic) 
interfacial width [104]. The chains in the bond fluctuation model comprise = 64 
segments which corresponds to an invariant degree of polymerization, JV = 240. 
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The time scale is set by r = R^/D = 1.5-10'' Monte Carlo steps, where each 
segment attempts on average one random, local displacement during a Monte Carlo 
step. In the simulations we use a cubic system of size L = 6.35i?e and average the 
results over 64 independent reahzations of the temperature quench. 

Fig. 13 compares the time evolution of the dynamic structure factor as calculated 
by dynamic SCF theory using a local Onsager coefficient (panel a) and non-local 
Onsager coefficient (panel b) with the results of the Monte Carlo simulation. As 
qualitatively expected from Cahn-HiUiard-Cook theory [101], there is a well-defined 
peak in the structure factor S'(q, t) 



that exponentially growth with time. For these calculations we have omitted the 
thermal noise. Therefore, composition fluctuations with q < Qc grow, while those 
with q> qcSie damped, and the structure factor for different times exhibit a common 
intersection at qc. No such well-defined critical wavevector qc is observed in Monte 
Carlo simulations or experiments, where thermal fluctuations limit the decay of 
fluctuations with q > qc- Clearly, there is a sizeable difference between dynamic 
SCF calculations using a local Onsager coefflcient and those which use a non-local 
A corresponding to Rouse-like dynamics. The peak of the structure factor in the 
dynamic SCF calculations with the local Onsager coefficient grows much too fast, 
and also the peak in the calculations occurs at larger wavevectors compared to the 
Monte Carlo simulations. Apparently, the calculations with the non-local Onsager 
coefficient agree much better with the simulation results. 

To quantify this observation, we investigate the growth (or decay) of density 
modes. The data in Fig. 14 show that the time evolution during the early stages of 
phase separation indeed is exponential and one can define a growth rate. The results 
for a local Onsager coefficient are presented in Fig. 12 and are bracketed by the 
approximations for the weak and strong segregation limits. 

In Fig. 15 we corroborate that the results of the dynamical SCF calculations using 
the non-local Onsager coefficient (114) agree with the EPD using the local Onsager 
coefficient (126) for the field W. As detailed in Sec. 5, however, the EPD calculations 
are computationally much less demanding. 

The growth rates for the dynamic SCF theory with local Onsager coefficient are 
displayed in Fig. 16 a) and compared to the results of the Monte Carlo simulations. 
The maximal growth rate occurs to too large wavevectors and the growth rate for 
q < qc is too large. The Cahn plot, R{q) /Dq^ vs. q^ presented in the inset, shows a 
linear behavior (cf. Eq. (132) in contrast to the simulations. 

Much of the discrepancy between calculations and Monte Carlo simulations is 
removed when we employ a non-local Onsager coefficient or EPD. For q < qcWe 
obtain almost quantitative agreement for the growth rate at early times. In the Monte 
Carlo simulations, however, deviations from the exponential growth is observed 
earUer than in the calculations and there is no decay of fluctuations for q > qc in 
the simulations. 




(134) 
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Abbildung 14. Several density modes displayed on a logarithmic scale versus time. The 
results were obtained through DSCFT calculations in a three dimensional system of length 
Lx — Ly — Lz = 6.357ie using 7x7x7 functions for a quench from = 0.314 to x-^ = 5. 
The expected exponential behavior during early stages of demixing is well reproduced. From 
Ref. [28]. 




Abbildung 15. Growth rates obtained through two dimensional DSCFT calculations using 
local and Rouse dynamics and EPD calculations. The DSCFT results using the pair-correlation 
function of a homogeneous melt and the EPD results are in good agreement. From Ref. [28]. 

Both effects can be captured if we include fluctuations in the EPD calculations. 
The results of EPD simulations including fluctuations are presented in Fig. 17. For 
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Abbildung 16. Corresponding relaxation rates to Fig. 13. Panel a) compares the Monte 
Carlo relaxation rates with DSCIT calculations with local dynamics. Panel b) compares 
Monte Carlo results with EPD and DSCFT calculations with Rouse dynamics. For earlier 
times good agreement in the growth region is found but Monte Carlo simulations show an 
earlier change in the exponential behaviour. From Ref. [28]. 




Abbildung 17. Relaxation rates obtained through Monte Carlo simulations and EPD 
calculations with random fluctuations for different time intervals. Fluctuations lead in both 
methods to an earlier change in the exponential behaviour of the increasing modes. From 
Ref. [28]. 

q < Qc and early times the calculations with and without fluctuations mutually agree, 
and resemble tiie simulation results. At later times, tiie growth rate decreases in the 
EPD and the simulations. For q > qc, fluctuations do not decay in the EPD but rather 
adopt their finite equilibrium strength. The growth rates extracted from the EPD for 
q > qc ratiier mirror the fact that some modes grow and other shrink during the 
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time window of the simulations, although the results have been averaged over 64 
independent, two-dimensional EPD calculations. 

The influence of the chain connectivity on the dynamics of the composition 
fluctuations does not only influence two-point correlation functions Uke the global 
structure factor, but it is also visible in the time evolution of composition profiles in 
the vicinity of a surface [105]. A description of the kinetics of phase separation is 
important because many multi-component polymer systems do not reach equilibrium 
on large length scales and the morphology of interfaces or self-assembled structures 
depend on processing conditions. 

The example highlights that both an accurate free energy functional as well as 
an Onsager coefficient that describes the dynamics of the molecules are important 
for a quantitative description. Of course, this study has concentrated on the simplest 
possible scenario. The assumption of the chain conformations being in equilibrium 
with the instantaneous composition field and the neglect of hydrodynamics have to 
be address to explore the behavior of more complex systems. 

8 

Summary and outlook 

We have briefly reviewed methods which extend the self-consistent mean-field 
theory in order to investigate the statics and dynamics of collective composition 
fluctuations in polymer blends. Within the standard model of the self-consistent 
field theory, the blend is described as an ensemble of Gaussian threads of extension 
Re. There are two types of interactions: zero-ranged repulsions between threads of 
different species with strength and an incompressibiUty constraint for the local 
density. 

In the EP theory one regards fluctuation of the composition but still invokes 
a mean-field approximation for fluctuations of the total density. This approach is 
accurate for dense mixtures of long molecules, because composition fluctuations 
decouple from fluctuations of the density. The energy per monomer due to composition 
fluctuations is typically on the order of x^bT ~ ksT/N, and it is therefore much 
smaller than the energy of repulsive interactions (on the order lA'sT) in the polymer 
fluid that give rise to the incompressibility constraint. If there were a coupling 
between composition and density fluctuations, a better description of the repulsive 
hard-core interactions in the compressible mixture would be required in the first 
place. 

Both the statics and the collective dynamics of composition fluctuations can 
be described by these methods, and one can expect these schemes to capture the 
essential features of fluctuation effects of the field theoretical model for dense 
polymer blends. The pronounced effects of composition fluctuations have been 
illustrated by studying the formation of a microemulsion [77]. Other situations where 
composition fluctuations are very important and where we expect that these methods 
can make straightforward contributions to our understanding are, e. g., critical points 
of the demixing in a polymer blend, where one observes a crossover from mean 
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field to Ising critical behavior [48, 49], or random copolymers, where a fluctuation- 
induced microemulsion is observed [62] instead of macrophase separation which is 
predicted by mean-field theory[61]. 

The apphcation of the dynamic SCF flieory [92] or EPD [30, 28, 105] to 
the collective dynamics of concentration fluctuations and the relation between the 
dynamics of collective concentration fluctuations and the single chain dynamics is 
an additional, practically important aspect. We have merely illustrated the simplest 
possible case - the early stages of spontaneous phase separation within a purely 
diffusive dynamics. In applications hydrodynamic effects [106, 107], shear and 
viscoelasticity [108] might become important. Even deceptively simple situations 
- like nucleation phenomena in binary polymer blends - still pose challenging 
questions [ 1 09] . Also the assumption of local equilibrium for the chain conformations, 
which allows us to use the SCF free energy functional, has to be questioned critically. 
Methods have been devised to incorporate some of these complications [73, 94, 108, 
107], but the development in this area is still in its early stages. 

Field theoretical simulations [71, 72, 77] avoid any saddle point approximation 
and provide a formally exact solution of the standard model of the self-consistent 
field theory. To this end one has to deal with a complex free energy functional as a 
function of the composition and density. This significantly increases the computational 
complexity. Moreover, for certain parameter regions, it is very difficult to obtain 
reliable results due to the sign problem that a complex weight imparts onto thermodynamical 
averages [77]. We have illustrated that for a dense binary blend the results of the 
field theoretical simulations and the EP theory agree quantitatively, i. e., density and 
composition fluctuation decouple [77]. In other systems, however, density fluctuations 
are of crucial importance. Let us mention two examples: (i) Field theoretical simulation 
methods can be used to study semi-dilute solutions [72, 74], i. e., a compressible 
ensemble of Gaussian threads with repulsive interactions. In dilute solutions the 
chains adopt a self-avoiding walk behavior, while in semi-dilute solutions the excluded 
volume interactions are screened on large length scale [16]. The very same model 
has been extensively studied by series expansions [110] and renormalization group 
theory [111] and those results provide an excellent testing bed to gauge the efficency 
of the field theoretical simulation method, (ii) It is a challenge to apply field 
theoretical simulations to particle-based models. This might provide valuable insight 
into the relation between coarse-grained parameters, Rg and xiV, of the standard 
model of the self-consistent field theory and the parameters (temperature, pressure/density 
and chain length) of a particle-based model. 
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